DFT-based transport calculations, Friedel's sum rule and the Kondo effect 
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Friedel's sum rule provides an explicit expression for a conductance functional, G[n], valid for 
the single impurity Anderson model at zero temperature. The functional is special because it does 
not depend on the interaction strength U. As a consequence, the Landauer conductance for the 
Kohn-Sham (KS) particles of density functional theory (DFT) coincides with the true conductance 
of the interacting system. The argument breaks down at temperatures above the Kondo scale, near 
integer filling, nda ~ 1/2 for spins a= \\.. Here, the true conductance is strongly suppressed by the 
Coulomb blockade, while the KS-conductance still indicates resonant transport. Conclusions of our 
analysis are corroborated by DFT studies with numerically exact exchange-correlation functionals 
reconstructed from calculations employing the density matrix renormalization group. 



The ground state density functional theory (DFT) 
owes its success to the fact that it proves enormously use- 
ful in the prediction of electronic properties of molecules, 
solids and surfaces. [2 [2] Therefore, applications towards 
the electronic transport properties of single molecules 
and self assembled monolayers came together quite nat- 
urally with the corresponding experimental successes in 
the field. [3J 13] Besides being useful for quantitative cal- 
culations, more fundamental properties of DFTs and the 
corresponding exact functionals have been an issue of in- 
tensive research. As an example we mention the ques- 
tion, what exactly is the nature of the approximations 
when using ground state DFT in combination with the 
Landauer formalism for transport calculations. E] 

The Landauer approach formulated in terms of non- 
equilibrium Green's functions [3 [8] and its validity for 
Kohn-Sham particles is also our topic in this work. It re- 
lates a single particle Hamiltonian, H K °, to the conduc- 
tance via the transmission function, G = T KS (£p)e 2 /h. 
The definition is 

T KS (E) = Tr T^ S (E)G(E)T^ S (E)G\E) (1) 

where r« s = i(E a - EjJ, G = (E - H KS - Z n - 
and the trace is over the Hilbert space associated with 
H KS . The self energies E Q describe the coupling of the 
KS-system to external reservoirs, a = £,1Z. They are 
given by a golden-rule expression, Ti a (E) — \V\ 2 g^ (E), 
where V is the coupling matrix element and g£° is a 
Green's function of the leads. [9] 

Often quantum transport is dominated by a single or- 
bital of the molecule or the quantum dot (QD), only. 
Therefore, in studies of correlated electron transport in- 
teracting level models are standard, e.g., the single im- 
purity Anderson model (SIAM, see also (|5|) 

H Q u=(- d N d + U (h dt - i) U di - i) , (2) 



where N d = h d ^ + h d ^ with h d lJ = d\d a and spin <r =f 
,4-. In such models an analogue version of is valid 
featuring retarded and advanced Green's functions of the 
interacting system (T a<acr > = T a 5 aa i), 

T{E) = -^^A d (E) (3) 

where we have introduced the spectral function of the 
interacting QD, A d (E)=iTr (7 (G (E)—G a (E)). [10] 

By comparing Eq. ^ and |T]) one might suspect, that 
in order to accurately reproduce the true value for the 
transmission, T(sp) = T KS (ep), it is neccessary for the 
KS-thcory to also reproduce the true spectral function, 
A d (E). It is easy to see, that the latter is not possible, 
however, unless U—0. To this end we recall that A d (E), 
carries the two Hubbard peaks at energies e d ± t//2.[TT] 
These peaks are not seen by the KS-system because the 
model ([2]) does not exhibit magnetism, so both spin 
channels, (t,4), are equivalent. Therefore, H KS is di- 
agonal, iJ^ s = e^ s S acr ', and the KS-spectral function, 
_4 KS = iTr cr (G — Gt), supports a single peak centered 
about £j , only. Despite of the absense of the Hubbard 
peaks in _4 KS we argue that in addition to the ground 
state density also the KS-conductance coincides with the 
true value of the interacting system, T = T KS . 

Moreover, we maintain that this statement is correct 
even though H KS is not unique in the sense that the 
exchange-correlation on-site potential, v^ s = e^ s — e d , 
can be complemented by an XC-contribution to the cou- 
plings, V —> V KS , as well. Different combinations 
(V ,v^ s ) produce an effective single particle Hamilto- 
nian with the correct ground state density. In fact, as will 
be demonstrated below, the on-site exchange-correlation 
(XC) potential and the coupling to the leads can be dras- 
tically different. Correspondingly, the resonance posi- 
tion, £j , and broadenings, Tc,ti, that determine ,A KS will 
be strongly XC-functional dependent. We will see that 
despite of this ambiguity in _ff KS the KS-conductance is 
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an observable that takes a unique value coinciding with 
the true conductance. 

We start with a general recollection about features 
of ground state DFT that follow directly from the first 
Hohenberg-Kohn theorem. [T^] According to this the- 
orem, we can reconstruct the external potential, f ex ( r ) 
that an interacting gas of N electrons is exposed to, if 
the ground state density n(r) together with the Hamilto- 
nian H in the absense of any u cx (r) is known. Thus we 
reproduce the full Hamiltonian H — Hq + V ex from our 
knowlege of n(r) (up to a constant shift in energy). 

Knowing H we can calculate, in principle, all equilib- 
rium response functions of the iV-particle system, pro- 
vided that the ground state is unique, once n(r) has 
been specified. Hence, we can consider such correla- 
tors to be a functional of the ground state density. [H] 
For example the dynamical density-density susceptibility, 
x[n](r, r', t — t'), can be thought of being such functional 
of n(r). Correlators can also be calculated at finite tem- 
perature. Again, they can be thought of as a unique 
functional of the ground state density under the condi- 
tion that the equilibrium state does not exhibit a broken 
symmetry which is not already encoded in n(r). 

A correlation function of special interest to us is the 
current-current correlation function, i.e. the conductiv- 
ity. Its longitudinal part describing the response to po- 
tential gradients is closely related to the density suscep- 
tibility: 



kc>x[ n ]( r 7 r', ui) 



d d 



-cr e [n}(r,r',uj) 



(4) 



dr dr 

The susceptibility, in turn, may be expressed as a conduc- 
tance in transport experiments that operate with elec- 
trodes for which the Fermi-liquid description holds true. 
Hence, in such geometries also the conductance is a func- 
tional of the ground state density: G[n]- Hence, this con- 
ductance can be calculated using the ground state DFT 
if the proper functional Q[n] is being used. 

The previous statement is as correct as it is useless for 
practical purposes unless a good approximation for Q[n] 
can be given. Of course, even if such an approximation 
would be known, in practice calculations would still suffer 
from inaccuracies in approximate XC-functionals used to 
obtain n(r). 

Friedel's sum rule for the SIAM. The complete 
definition of the SIAM-Hamilonian reads [13] 



H — Hqd + ^ Ha + Ht 

a=C,TZ 

M-X 
x—l,a 



(5) 



H T = -V 



<■■ ) ■ (6) 

(7) 



The Cx}a,a denote fermionic annihilation (creation) oper- 
ators at site x, lead a = 1Z, C The model is of interest to 



us because it affords the Friedel sum rule|13|. 

sin 2 [7rn dcr (£ F )] 



^dcr(eF) 



T/2 



(8) 



where T = Tc + T^. By combination with ([3]) we relate 
the particle number on the quantum dot, = n^f = 

ndi, to the conductance, 
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2e 2 TrT 



h (r/2)^ 



sin 
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(9) 



The identity constitutes an exact analytical expression 
for a conductance functional C?[n].It is remarkable that 
in the case of symmetric coupling, Tc—I^tz, Eq. ([9| re- 
lates the conductance associated with H to a single sys- 
tem characteristics, only, which is the particle number 
on the QD, N^. This implies, that any change in the pa- 
rameters of H leaves the conductance invariant, provided 
that Nd is unchanged and that the conditions of appli- 
cability of Q are still valid. One requirement for this 
is that both leads are (effectively) non-interacting and 
free of backscattering, so that the Fermi-liquid picture 
holds a zero temperature. The last condition is impo- 
rant because else ir-rido- would not properly account for 
the number of bound states (per spin) and the scattering 
phase. 

Friedel's sum rule and ground state DFT. The KS- 
Hamiltonian of ground state DFT reads 



H 



KS 



Na 



E 



V v KS 



N,. 



H T (10) 



where the onsite XC-potentials, eJp[A^] = e d + w d s [iV] 
and v^JyN], axe functional of the local particle density 

N x . At first sight, the KS-Hamiltonian, H KS , does not 
meet the requirements for the validity of (|9|, since the 
XC-functional extends into the leads and thus might con- 
tribute to the scattering phase shift: 5(e-p) — > (5 KS (ep)- 
However, this impression is misleading as can be seen 
from the following argument. The Friedel sum rule in 
its general form relates the extra scattering phase shift, 
S(e-p), induced by decreasing £d (down from oo) to the 
extra spectral weight, AA(E), thus generated: 



S(s F ) 



dE 
2tt 



(11) 



The right hand side denotes the total change in the parti- 
cle number per spin, A N a = J2 X AN X(7 , associated with 
the occupation of this extra weight. A special aspect of 
Eq. Q is the implication that the sum rule is exhausted 
already by the density change on the QD, AN a = 
the net charge accumulated by bringing £d down is en- 
tirely concentrated in the QD. This latter statement must 
be correct also for the KS-theory, by definition. Since 



also KS-particles obey the relation (11) we conclude: 
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FIG. 1: Ground state density per spin of a QD (at 2 = 0) 
coupled to a non-interacting reservoir with M=100 sites for 
growing on-QD interaction U = 0.0, 0.6, 4.0. Parameters: 
£d=0.2, V = 0.3, bandwidth of conduction electrons: 2t = 2. 



FIG. 2: XC-correlation potential corresponding to the evolu- 
tion of the density shown in Fig. [T] The on-dot potential is 
denoted v$ s = ef s . 



5(e F ) = S KS (e F ) = 7riV d /2 and Eq. ^ holds for them as 
well despite of the presence of u^f in the leads. 

Functional ^ is valid at any value of the interaction 
strength U and at temperatures below the Kondo scale, 
T K ~ 2t^2T/Ue~ u/ir near integer filling N d « l.[13] In 
this context the Abrikosov-Suhl resonance underlying the 
Kondo-cffect plays a crucial role. At temperatures above 
Tk it is not developed and the conductance is strongly 
suppressed due to the Coulomb blockade; relation (J9j) is 
strongly violated. The XC-functional of DFT must be 
very sensitive to Kondo-physics. This is obvious for the 
following reason: we have seen that KS-transport repro- 
duces the exact transmission. To reproduce the resonant 
transport (Kondo) scenario in the regime N d ~ 1, the 
KS-level of the quantum dot must be half filled for each 
spin, implying \e F — £^ s | < T in the Kondo regime even 
if the bare position of the level — e d 3> T. 

To illustrate and extend our analysis we have calcu- 
lated the ground state density and corresponding exact 
XC-functionals employing the density matrix renormal- 
ization group method and "backward" DFT. The ap- 
proach has proven useful before in the context of the 
interacting resonant level model (IRLM). [2] We have 
adapted our technology here to treat the SIAM. [T3] . 

In our calculations we consider a coupling to a single 
lead only in order to reduce the computational effort. As 
long as the ground state is concerned the case with two 
leads and symmetric couplings, Tc = T^, has an exact 
mapping into the single lead case, essentially because the 
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odd combination of tunneling operators, 

c£\)/\/2 decouples from the QD. 

In Fig. [I] we display the evolution of the density in 
a QD with e d slightly above e F . Without interactions, 
the QD is empty. Upon increasing U, the dot fills be- 
cause in the spirit of a Jellium model we have defined 



the interaction in ^ with respect to density fluctuations 
against a background n^ a = 1/2. The density in the 
leads, x > 1, exhibits the typical Friedel-oscillations with 
their 2k F periodicity and the l/x-envelop. Their ampli- 
tude is controlled by the boundary condition which is set 
by the QD. Since its occupation changes from Ad ~ 
to N d w 1 the phase shift associated with backscattering 
increases by 2-nn da . This is why at large U the Friedel- 
oscillations are anti-phase with the case U =0, inset Fig. 
0. We witness a signature of Kondo-physics, cf. [IB] . 

Fig. [2] shows the KS-potential on the QD, e d cs and in 
the lead, v£ s , corresponding to the evolution of the den- 
sity, Fig. [IJ On the dot v^ s shows the expected behavior 
from e^P = e d at U — to e^ s ~ at large interaction, 
here U = 4. In the leads the potential oscillations fol- 
low the Friedel-oscillations of the density and introduce 
the interaction corrections. The oscillation amplitude de- 
pends in a non-monotonic way on U, increasing from the 
non-interacting fixpoint and decreasing again when ap- 
proaching the strong-coupling, Kondo fixpoint, see inset. 

In a KS-theory based on a local KS-potential, the hy- 
bridzation V coincides with the one from the original 
coupling H T . This suggests that the on-dot spectral func- 
tion, A KS (E) of the KS-system has a width close to the 
non-interacting one, V w T KS . Fig. [3] fully supports this 
point of view. In addition, it also shows that indeed there 
is only a single maximum (Hubbard peaks don't exist in 
KS-theory) which will be getting closer to e F with in- 
creasing interaction strength U. 

Fig. [4] compares exact transmission curves obtained 
via the Friedel sum rule Landauer's formula (Tc = T^) 



T(e F ) 



(172) 5 



£F )2 + (T/2Y 



(12) 



The good agreement between the results obtained with 
both methods illustrates the point emphasized above: 
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FIG. 3: Local spectral function of the KS-system, A d = 
27rpi mp , on the QD for interactions U = 0,2. Parameters: 
e d = 0.1, V = 0.15..., £ F = 0. 



Even though the KS-spectral function is not physical, 
the associated transmission is close to the precise value. 

We briefly mention that this point can be highlighted 
further by constructing a version of KS-theory, in which 
not only the diagonal elements of the density matrix, i.e. 
n xa but also off-diagonal elements are faithfully repro- 
duced. This can be achieved by adding to the on-site po- 
tential also a modification of hopping matrix elements, 
i.e. V — > V KS such that the expectation value of the 
kinetic energy of the impurity coupled to the first lead 
site within the KS description matches the one obtained 
from DMRG; for the technical details see [15]. Respond- 
ing to this change the width of the KS-spectral function 
of the modified theory is no longer close to the original 
one, r KS — > r KS . Since the charge in the QD must re- 
main unchanged, we expect a compensating shift in the 



on-site energy, e 



KS 



-KS' 



As can be seen from Fig. 



[4] despite a substantial change, V = 0.3 is replaced by 
V = 0.1..., the transmission when evaluated via (12 1 is 
not changed. In view of Eq. (12 1 this finding is eas- 



ily understood: T(ep) is determined by the same ratio, 
r KS /(£j s — Sp), that also fixes the density. Since by con- 
struction in all KS-models the density is the same, we 
have r KS /( £ KS _ ep ) = r Ks' /(£ Ks' _ £p)j so the trans _ 

mission remains the same, also. 

Conclusions. We briefly discuss two generalizations 
of the proceeding analysis. First, our argument was as- 
suming a single level, only, while a real QD, e.g. a 
molecule, exhibits in general several levels. One expects, 
however, that the main conclusion remains valid as long 
as all levels contribute independently to the transport 
current. In particular, the sum rule ^ should remain 
a useful approximation for the true density functional. 
Second, our analysis heavily relies on the Kondo effect 
restoring full transmission in the case of single occupa- 
tion of the dot. Suppose, that the Kondo temperature is 
very low, and that the measurement is done at slightly 
higher temperatures. Then transport is dominated by 
the Coulomb blockade and the conductance is strongly 
suppressed, T(ef) ~ (r/C/) 2 1. On the other hand 
the particle density is essentially still the ground state 
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FIG. 4: Comparison of conductances calculated via the 
Friedel sum rule (FSR) and directly from the Landauer for- 
mula for KS-particles, Q. Parameters V = 0.3, U = 1,8,2.1 
are considered. 



one, n(r), i.e. it is largely insensitive to this change and 
in particular ps 1/2. We conclude that in this case 
the Friedel sum rule ^ does not hold and that KS-theory 
(without breaking spin rotational invariance) does not re- 
flect this change, i.e. we still have: T KS (ep) ~ 1. 
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